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ABSTRACT 



Four different light signals, all of the same optical frequency, 
but, having different states of polarization, are transmitted 
through an optical device and the mean signal delay of each 
of the light signals is measured. Calculations, based upon the 

relationship, x g **x 0 - l /i Q s, describing the polarization 
dependence of x g (a measured mean signal delay) through 
the device as a function of x 0 (a polarization independent 

delay component of the device), If (the PMD vector at the 
device input) and s (the input Stokes vector of the light 
signal), yield the PMD of the device. Also, by comparing 
data taken at adjacent wavelengths, the chromatic dispersion 
of the optical device can be accurately measured even in the 
presence of PMD. 

15 Claims, 3 Drawing Sheets 




LAUNCH 
POLARIZERS 



04/13/2004, EAST Version: 1.4.1 



U.S. Patent Feb. 11,2003 Sheet 1 of 3 US 6,519,027 Bl 



FIG. 1 



MODULATED 




POLARIZATION 


SOURCE 




CONTROLLER 




FIG. 2 




LAUNCH 
POLARIZERS 



04/13/2004, EAST Version: 1.4.1 



U.S. Patent Feb. ll, 2003 Sheet 2 of 3 



US 6,519,027 Bl 



80- 



FIG. 3 



to 

Ql 



3 



60- 



< 40- 

o 
to 



i 

LxJ 



20- 



INPUT POLARIZATION: 



€ Si 

a 60° LINEAR 
* 120° LINEAR 
S3 




1543.75 



1544.0 



WAVELENGTH (nm) 



1544.25 



04/13/2004, EAST Version: 1.4.1 



U.S. Patent Feb. 11,2003 Sheet 3 of 3 US 6,519,027 



FIG. 4 




1543.8 1544.0 1544.2 



WAVELENGTH (nm) 



04/13/2004, EAST Version: 1.4.1 



US 6,51 

1 

POLARIZATION MODE DISPERSION 
MEASUREMENT 

BACKGROUND OF THE INVENTION 

This invention relates to optical measurements, and par- 
ticularly to the measurement of the polarization mode dis- 
persion of optical devices, particularly optical fibers used in 
telecommunication networks. Although the polarization 
mode dispersion measurement can be made at a single light 
wavelength, by making the measurements at plural 
wavelengths, other optical parameters, e.g., chromatic 
dispersion, can be determined. 

Polarization mode dispersion (PMD) is a distortion 
mechanism (like chromatic dispersion) that causes optical 
devices, such as single-mode fibers, optical switches and 
optical isolators, to distort transmitted light signals. The 
relative severity of PMD (which is a function of the wave- 
length of the transmitted light) has increased as techniques 
for dealing with chromatic dispersion have improved, trans- 
mission distances have increased, and bit rates have 
increased. Negative effects of PMD appear as random signal 
fading, increased composite second order distortion and 
increased error rates. 

PMD is due to differential group delay caused by geo- 
metrical irregularities and other sources of birefringence in 
the transmission path of the optical device. For example, a 
single-mode fiber (SMF) is ideally a homogeneous medium 
supporting only one mode. In practice, it supports two 
propagation modes with orthogonal polarizations. When a 
lightwave source transmits a pulse into a SMF cable, the 
pulse energy is resolved onto the principal states of polar- 
ization of the fiber. The two groups of pulse energy propa- 
gate at different velocities and arrive at different times 
causing pulse broadening and signal distortion. 

Determining the PMD of in-place (installed) optical fibers 
is useful for determining the capacity of the fibers for 
transmitting new telecommunication services of ever 
increasing bandwidth requirements and for the design and 
control of PMD compensators. Also, because the PMD of 
any given fiber can change rapidly with time (owing to 
random asymmetrical stresses along the length of the fiber), 
repeated PMD measurements of the fiber can be used in a 
feedback system for compensating for variable PMD dis- 
tortions. 

The PMD of a fiber is commonly characterized by two 
specific orthogonal states of polarization called the principal 
states of polarization (PSPs) and the differential group delay 
(DGD) between them. This can be described at an optical 
angular frequency, co, by the 3-component Stokes vector, 

Q =At q , where q is a unit Stokes vector pointing in the 
direction of the faster PSP, and the magnitude, Ax, is the 
DGD. Typical DGD values encountered in transmission 
systems range between 1 (picosecond) ps and 100 ps. 

Known methods for determining PMD vectors include the 
Jones Matrix Eigenanalysis (JME) technique and the Miiller 
Matrix Method (MMM). Each of these techniques uses a 
tunable, continuous-wave laser and a polarimeter to measure 
the output polarization states for two (or three) different 
input polarization launches at two optical frequencies. The 
PMD vector is then calculated for the midpoint frequency. In 
addition to determining the output PMD vector, the Miiller 
Matrix Method determines the rotation matrix of the fiber at 
each frequency and thus the input PMD vector can be 
calculated. 
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Shortcomings of these techniques are that they are some- 
what difficult to implement, particularly in a feedback 
system, because they require frequency differentiation of 
measured data at plural optical frequencies. As described 

5 hereinafter, an advantage of the method according to the 
present invention is that it requires measurements of data at 
only one optical frequency. Also, while, as hereinafter 
described, measurements can be made at plural optical 
frequencies, separate, single frequency measurements are 

10 made that do not involve frequency differentiation. 

A further advantage of the invention is that it can be used 
to accurately measure chromatic dispersion of an optical 
device in the presence of PMD. Other known methods for 
measuring chromatic dispersion tend to give inaccurate 

1 5 results if the devices have PMD. 

SUMMARY OF THE INVENTION 

At the input of an optical device, typically an optical fiber 
2Q link within an optical fiber telecommunications network, 
four different light signals, all of the same optical frequency, 
but having different states of polarization, are transmitted 
along the fiber and the mean signal delay of each of the fight 
signals is measured. By calculation (described hereinafter), 
1S the four mean signal delay measurements yield the desired 
PMD vector of the optical device. Additionally, by repeating 
the mean signal delay measurement at multiple optical 
frequencies (i.e., at a different optical frequency for each set 
of time delay measurements) determination can be made of 
30 the first and higher-order intrinsic (polarization-averaged) 
chromatic dispersion of the device being measured, as well 
as the higher-order PMD of the device. 

DESCRIPTION OF THE DRAWING 

35 FIG. 1 is a schematic illustration of a system for imple- 
menting the present invention; 

FIG. 2 is similar to FIG. 1 but showing a more specific 
example of a suitable system for implementing the inven- 
tion; 

40 FIG. 3 is a graph showing the wavelength dependence of 
the signal delay of each of four fixed input polarization 
signals transmitted through a single mode fiber; and 

FIG. 4 is a graph showing the wavelength dependence of 
the input PMD vector of the fiber obtained using the data of 
FIG. 3. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

50 As above-described, PMD is usually (and herein) 
described as a three dimensional vector having a magnitude 
and direction, i.e., having three components. Both an input 
and an output PMD vector can be determined for any optical 
device. The PMD vectors are a function of the wavelength 

55 of the light transmitted through the device. The present 
invention involves a new method for making measurements 
of an optical device and determining the three components 
of the various PMD vectors. The measurement and calcu- 
lations made are based upon the known (Mollenauer, L. F. 

60 and Gordon, J. P., Optics Lett., Vol. 19, pp. 375-377, 1994) 
relationship, 

T g *r 0 -YiQs (1) 

65 where s and Q are, respectively, the Stokes vectors of the 
light and the PMD vector at the fiber input. Equation (1) 
describes the polarization dependence of the mean signal 
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delay, x^ through the fiber as defined by the first moment of 
the pulse envelope in the time domain with x 0 being a 
polarization independent delay component. Eq. (1) assumes 

that x 0 and * (co) do not vary significantly over the 
bandwidth of the signal. 

The definition of mean signal delays involving 
"moments" of the output signals is known and described, for 
example, in the afore -cited publication by Mollenauer and 
Gordon as well as by Elbers, J. P., et al, "Modeling of 
Polarization Mode Dispersion in Single Mode Fibers/ 1 
Electr. Lett., Vol. 33, pp. 1894-1895, October 97; Shieh, W., 
"Principal States of Polarization for an Optical Pulse/' IEEE 
Photon. Technol. Lett; Vol. 11, No. 6, p. 677, June 99; and 
Karlsson, M., "Polarization Mode Dispersion-Induced Pulse 
Broadening in Optical Fibers," Opt. 

Lett., Vol. 23, pp. 688-690, *98. 

More precisely, mean signal delay, x g , is expressed as the 
difference of the normalized first moments at fiber output 
and input, 



. wiU)-yyi(Q) 



Q 2 -2(i 



■v) 



10 



20 



(2) 



where z is the distance of propagation in the fiber. Here 

W-Jdt Ef E-JdcoEfE is the energy of the signal pulse 

represented by the complex field vector E (z,t) with Fourier 

transform E(z,co), and W 1 (z)=/dttEf E =j JdcoEfE^ is the 
first moment. 

Because, in equation (1), there are four unknowns, i.e., x 0 
(fiber intriasic group delay) and the three components of the 

PMD input vector, Q, four measurements are necessary to 

determine Q and at a given wavelength. Each set of 
measurements comprises launching a set of four light signals 
of respectively different polarization states into the fiber and 
measuring the mean signal delay of each of the signals by 
known means. 

A variety of different fight polarization states can be used, 
but for simplicity of the mathematical analysis, a convenient 
set of polarization states are those coinciding with the 
Poincare* sphere axes, S^-nS^ S 2 , and S 3 , i.e., three different 
linear Polarization states and a circular polarization state. (A 
"Poincare* sphere" being a graphical representation of all 
possible Polarization states on a surface of a sphere where 
each point on the sphere represents a different polarization 
form; see, for example, W. Shurcliff, Polarized Light: Pro- 
duction and Use, p. 16, Harvard University Press, 
Cambridge, Mass., 1962, A specific example of selected 
polarization states is provided hereinafter.) 

By identifying the measured mean signal delays for each 
polarized state as x gl9 t^, and x^ respectively, the 
calculated components from equation (1) are as follows: 

to-^(^i«,(-i)) (3) 

and the three components of the fiber input PMD Stokes 
vector Q are: 



35 



polarizations, can be generalized to any four non-degenerate 
input polarization states.) 

As above-described, implementation of the present inven- 
tion involves the launching of a set of different test light 
signals into a device under test, e.g., an optical fiber, and 
measuring the mean signal delay for each signal. 

The basic technique according to the invention is illus- 
trated in FIG. 1, which shows a test system including a 
source 10 of modulated light at a single optical frequency 
and a polarization controller 12 for adjusting the polarization 
of the light from the source 10 to four successively different 
polarization states which are successively introduced into 
and transmitted through an optical device 14 under test, e.g., 
a length of optical fiber wound on a coil. The light signals 
exiting the test device 14 are successively fed into an 
apparatus 16 for measuring the mean signal delay for each 
test signal transmitted through the device. The measured 
data provides, by equations (3) and (4), the polarization 
independent delay and the input PMD Stokes vector of the 
device under test at a single optical frequency. The mea- 
surements are repeated at successively different optical 
frequencies for determining the frequency dependence of 
these factors. Also, the data taken at multiple frequencies 
allows the determination of the chromatic dispersion of the 
device under test. This is described hereinafter. 

As noted, the system shown in FIG. 1 measures the PMD 
vector at the input to the optical device 14. An advantage of 
the system shown in FIG. 1, in comparison with the prior 
mentioned known Jones (JME) and MuHer (MMM) tech- 
niques for measuring input PMD vectors, is that no output 
polarimeter is required. Conversely, if the output from the 
optical device is passed into a polarimeter, the rotation 
matrix of the optical device, R(co), can be determined. Then 

the PMD vector at the output of the optical device, Q^w), 
can be determined from 



50 



55 



(4) 



(Described hereinafter is how the above mathematical 
relationships, for the "special" Poincare sphere axial 



R,(f>)=.K(to)R(<i>), 



(5) 



where Q(co) is the input PMD vector from eq. (4). For 
instance, for the previously mentioned launches coinciding 
§ 2 , and S 3 , the matrix R(to) is given by 



45 



(6) 



where T 1; -, T^, and T 3j - (j»l, 2, 3) are the jth components of 
the output Stokes vectors measured by the polarimeter for 
launches S 19 § 2 , and § 3 , respectively. 

Apparatus for performing each of the functions indicated 
in FIG. 1 are generally known, and different specific 
arrangements for practicing the invention can be devised by 
persons of skill. FIG. 2 shows but one example of an 
arrangement useful for implementing the invention. 

Light of a desired wavelength is provided by a wave- 
length tunable laser 20 (Hewlett-Packard 81682A). A direc- 
tional coupler 22 (Etek SWBC2201PL213) splits the light 
from the laser into two parts, one of which enters a waveme- 
ter 24 (Hewlett-Packard 86120B) that measures the wave- 
length of the light, the other of which is coupled to an optical 
modulator (MOD) 26 (Lucent Technologies X2623C). An 
electrical signal from an oscillator 28 (1 GHz) (part of 
Hewlett-Packard 83420A) is also connected to the modula- 
tor. The optical modulator imposes a sinusoidal modulation 
on the optical power of the tight. This modulated signal is 
amplified by an optical amplifier (A) 30 (Lucent Technolo- 
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gies W1724CDDAD), passed through a tunable optical filter At»0.45 at the maximum DGD, where here co m is 2rc times 

32 (JDS Fitel TB 1500B) tuned to the wavelength of the the modulation frequency. FIG. 4 also shows data for the 

signal and enters a linear polarizer 34 (Etek determination of the chromatic dispersion of the test fiber 44 

FPPD2171LUC02). A polarization controller (PCI) 36 (FIG. 2). Chromatic dispersion is the change in the intrinsic 

(Fiber Control Industries FPC-1) is used to adjust the s (polarization-averaged) group delay x 0 through an optical 

polarization of the signal so that it nearly matches the j . ... , „ . . r f 

r ! . - e , . r% * t*u i ■ ia device with wavelength. For instance, if measurements are 

polarization of the linear polarizer 34. 1 nis polarizer 34, , , , * , * . . . * . 

which is not essential to the measurement, serves to make taken at two wavelengths, X x and X 2 , the chromatic disper- 

the results of the measurement insensitive to polarization slon D at thc average of these wavelengths, can 

changes occurring in the components and in the fiber pre- be obtained by dividing the change in x 0 by the change in 

ceding the polarizer. After the polarizer 34, the signal passes wavelength: 
though a launch polarizer box 40. This box 40, which 

functions by switching polarizers into or out of the signal i m - T o( A 2)- r o(*i) 0) 

path, selectively modifies the signal polarization from its ™ 1 * 2 " * " \ 2 -h 
input polarization to one of four output polarization states. 

Three of those states are provided by putting an appro- 15 

priate polarizer in the beam path while the fourth state is In FIG. 4, the chromatic dispersion of the fiber corresponds 

provided by switching all of the polarizers out of the beam to the slope of the curve plotting x 0 versus wavelength. For 

path and using the input signal polarization unmodified. This example, at a wavelength of 1544.0, the dispersion is 73.3 

input polarization can be adjusted to a desired state using n S ' nm 

polarization controller (PC2) 42 (Fiber Control Industries 20 ' _^ 

FPC-1). Alternatively, the fourth polarization could be pro- Also shown in FIG. 4 is a measurement of the input Q (co) 

vided by a fourth switchable polarizer in launch polarizer on t he same fiber using the aforementioned Miiller Matrix 

box 40. After passing through the launch polarizer box 40, Method for determining PMD vectors. The data for the M 

the signal is launched into the device under test, shown here » Uer method wag takeQ immediately before the measure . 

as a fiber span 44^ After ^ the fiber span the signal, now 25 ^ accordin t0 the t Nation. The two techniques 

modified by the PMD of the fiber span, enters a receiver 46 ^ maximum ^ of 7 

(part of Hewlett-Packard 83420A) which generates an elec- ^ v j v 

trical signal with a voltage proportional to the power in the While it is convenient to launch the specific polarization 

optical signal. A network analyzer 48 (Hewlett-Packard states, S^-Sj, S 2 , and S 3 , as previously noted, equations (3) 

8753D) measures the difference between the phase of the 30 and (4) can be generalized to any four non-degenerate input 

signal provided by the oscillator 28 and the phase of the polarization launches % (i«l,2,3,4) that span Stokes space, 

amplitude modulation of the received signal. Prior to enter- The corresponding signal delay measurements, x gi , can be 

ing the receiver 46, the signal can optionally be tapped using grouped into a four -dimensional (4-D) vector 

a directional coupler 50 (Etek SWBC2210PL213) and the x^x^x^x^x^), and similarly, the desired delay compo- 

tapped light sent to a polarimeter 52 (Hewlett-Packard 35 nen ts can be grouped into a 4-D vector x s =(x 0 ,Q 1 ,Q 2 ,Q 3 ) 

8509B). This polarimeter 52, which is not needed for the The four launches together with Eq. (1) set up the required 

inventive measurement techique, is used in the measurement f 01ir linear equations for the four unknowns represented by 

of PMD by a prior technique such as the Miller Matrix x, yielding x g =P x, where P is a 4x4 matrix he components 

Method. This allows comparison of the results obtained by 0 f me launched polarizations, 

the different techniques. The polarimeter 52 is also conve- 4 ° 

nient when setting the polarization controller 42 to a desired 2 -s u ~s n -s 13 (8) 
state. 

An example of measurements made according to the 
invention using the apparatus shown in FIG. 2 follows. 

Four polarization states (linear polarization states at 0°, 45 
60° and 120°, and a circular polarization state) were 
launched into a 52-km optical fiber span having a known 

average PMD of 35 ps and a chromatic dispersion of Here S if is the jth component of polarization launch S The 

approximately -6 ps/nm at 1555 nm. equivalent of Eqs. (3) and (4) is then-c=P _1 %g, and so t 0 and 

FIG. 3 shows the signal delays (vertical axis) for each 50 ^ co CQts of £ ^ be dctermmcdby iaver ting 

input polarization as a function of light wavelength ^ matrix p at each fr of intefest 

(horizontal axis). (Because the intrinsic signal delays /rx 

through the 52-ta fiber spaa are about 260 /*, "relative" A For this gf nera = ase ' ® can ^ be ^ l .° 

signal delays are plotted in HG. 3 where a constant delay has determme the output PMD vector when a polanmeter is 

been subtracted from the actual measured delays.) The PMD 55 inserted at the output In order to determ.ne the rotation 

of the fiber span causes the polarization dependence of the matrix of the flbe j; ^ an y two of "» four non^egenerate 

signal delay at each wavelength, and the variations in those input launches S ; can be used. We use S a and S b to 

relative signal delay curves ^originate from the change of designile ^ ^ chosen ^ ^ and as 

direction and magnitude of Q (co) as well as from under- 6Q the corresponding measured output Stokes vectors. We first 

lying chromatic dispersion. create a set of three orthogonal vectors at the input, § 2 '» 

FIG. 4 shows the input PMD vector Q (co) and the relative and S 3 \ where S 1 , -"? fl , ^«~S a xS b l\S a x~S b \, and § 2 '=§ 3 'x 

263 polarization -independent delays, x 0 , for this fiber, deter- §j\ Similarly we create a set of three orthogonal vectors at 

mined from the data of HO 3 Here, t 0 is the polarization- ou f f and f where f f f^fyf^ 

mdependent delay through the fiber span at each wavelength 65 ^ * 

compared to the delay at the first wavelength, 1543.70 nm. T a x T J, and ly«lyxtY. The rotation matrix of the fiber is 

The modulation frequency was 1 GHz, providing co m then R-R r Ry, where 



2 -5 U -S n -5 i3 

2 -S 2 i -S22 -Sis 

2 -S31 -£32 -£33 

2 S41 -542 -543 
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and 



^11 S'll ^13 

$11 S 'tt 

^31 ^32 5*3 



(9) 



(10) 



A first trial value for x 0 can be obtained by linearizing Eq. 
(15), 



Here S l7 -' 7 S 2j -', and S 3 y (j=l, 2, 3) are respectively the jth 
components of the set of three orthogonal vectors at the 
input, S x \ S 2 ', and S 3 \ while T iy ', T^, and T 3 / (j=l, 2, 3) are 
respectively the jth components of the set of three orthogo- 
nal vectors at the output, ly, ly, and ly. 

Instead of determining moments from pulse -shape 
measurements, the signal delays, % g£l can also be obtained (in 
approximation, i.e., to an acceptable degree of accuracy, 
e.g., 6% as described hereinafter) by observing phase shifts 
of the modulation of a sinusoidal amplitude modulated 
signal and benefiting from the precision of sensitive phase- 
detection techniques. (Note that Williams [Electron. Lett. 
Vol. 35, pp. 1578-1579, 1999] has used sinusoidal modu- 
lation for the determination of the scalar DGD.) For most 
purposes, it suffices to only measure changes in i gi with 
polarization and optical frequency and to not resolve the 
ambiguity presented by the use of a signal with periodic 
modulation. For a sinusoidal intensity -modulated signal 
with the assumption of frequency -independent PSP's and 
DGD, Eq. (1) becomes 

tan <i>„(^-t 0 )=- q - s tan((0,„AV2), (11) 

where co m is the angular modulation frequency. When the 
four launch polarizations coincide with the Poincare' sphere 
axes, §!,-§!, §2> anc * §3? Eq. (3) is still valid for sinusoidal 
modulation :*t 0 =! ^('V +T s(- ]))• Using the relation q 1 1 +q 2 2 + 
q 3 2 -l, expressions for the magnitude, At, and components, 

q,., of the unit vector, q , pointing in the direction of the 

PMD vector, Q , can also be determined and solved exactly. 

In the more general case we can use any four non- 
degenerate input polarization launches s , a, b, c ) that 
span Stokes space, where s 1 =S 1 , a 0 =(a 1? a^, a3), s fc =( D i> 

b 2 , b 3 ), and s C =(C U C 2 , C 3 ), and measure the correspond- 
ing delays x gl , x ga , x Mb , and x . We first express & 19 S 2 , and 

§ 3 , in terms of ~S a , !> b , and ~$ ct 

SrafT.+p/Tfc-KfcT^i-W), (12) 
where the coefficients a-, p„ and are obtained from 





A 


7i 




ai 


a 2 


ay 




02 


72 






b 2 


by 


ay 


Jh 


73 




C] 




Cy 



25 



30 



45 



50 



(13) 



55 



Substituting § 1 =o. i ~s a +$ 1 ~s b +y 1 ls c into Eq. (11) gives 

tan a> w (T x i-T 0 )— q -(aj s s 6 +y t s ^\&n(<* m to/2) y (14) 

leading to a transcendental equation for x 0 : 

tano) m (t, 1 -t 0 )=a 1 tancD m (t ta -T 0 +p 1 tana) M (T; i6 -i; Q )+ 

Tfitano^^,^. (15) 



60 



65 



(16) 



10 



After solving forx 0 and using qj 2 +q2 2 +q3 2 -l, an expression 
for the magnitude, Ax, can be found: 

tan 2 [w ffl Ar/2) = 0?) 
^ [Qf ; tantJ w (T^ - r 0 ) + p;tenu> m {T gb - r 0 ) +■ y.tan^fr^ - t 0 )] 2 . 



The components, q,-, of the PMD vector, Q, can then be 
determined from 



20 



gjtam^Ty, - t q ) + p;mti> m (r gb - r 0 ) 4- y.tanayfT^ - t 0 ) (18) 



Although the above procedure will yield any computational 
accuracy desired, it is often not necessary. For small modu- 
lation frequencies, co m , we can approximate tan (x)sx in Eq. 
HI), reducing to the earlier expressions (i.e. Eqs. (3) and 
(8)). These linear expressions are valid for sinusoidal modu- 
lation to within 6% as long as o) m Ax<^4. For instance, for 
the peak DGD we observed here, 70 ps, 6% accuracy will be 
obtained for modulation frequencies, fm=tom/2jt, less than 
1 .8 GHz by using the linear expressions. 

As above-explained, the inventive method involves mul- 
tiple launchings of different polarized light states, determin- 
ing the mean signal delays of the different polarization 
states, and calculating the PMD vectors. Basic techniques 
for measuring the group delay of an optical device are 
known, see, e.g., Y. Horiuchi, et al., "Chromatic Dispersion 
Measurements of 4564 km Optical Amplifier Repeater 
System," Electronics Letters, Vol. 29, pp. 4-5, Jan. 7, 1993. 
Additionally, by comparing data taken at adjacent wave- 
lengths with the inventive method, chromatic dispersion of 
optical devices can be accurately measured even with the 
presence of PMD. 

What is claimed is: 

1. A method of measuring the polarization mode disper- 
sion (PMD) at a first optical frequency of an optical device 
comprising launching at an input of said device at least four 
test light signals at said first frequency but at respectively 
different polarization states that span Stokes space, 
measuring, at an output of said device, the mean signal delay 
of each of said test signals, and calculating a three- 
component Stokes vector describing the PMD of said device 
at said input. 

2. A method according to claim 1 wherein the calculation 
is based upon the formula: 

T^o-HQ-S _ (1) 

where s is the input Stokes vector of a light signal trans- 
mitted through the device, Q is the input PMD vector for the 
device at the optical frequency of the light signal, x g is the 
mean signal delay, and x 0 is the intrinsic group delay of the 
device. 

3. A method according to claim 2 including determining, 
at said device output, the polarization state of each test light 
signal for determining the PMD of said device at said output 
thereof. 

4. A method according to claim 3 including using at least 
two of said test signals of known polarization state for 
determining the rotation matrix of said device. 

5. A method according to claim 4 wherein the measuring 
of the mean signal delays is accomplished in approximation 
by measuring phase delays of sinusoidally modulated sig- 
nals. 
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6. A method according to claim 3 wherein the measuring 
of the mean signal delays is accomplished in approximation 
by measuring phase delays of sinusoidally modulated sig- 
nals. 

7. A method according to claim 2 wherein said four 
polarization states are those coinciding with the axes 

S 2 , and § 3 of a Poincare* sphere representation of all possible 
light polarization states and wherein the calculated compo- 
nents of formula (1) are 

and the three components of said Stokes vector are: 

a 2 »2(T 0 -x i2 ) 
523-2(^-^3) 

where x gl7 T^^, arj d x g3 are tne respective measured 
mean signal delays for said test light signals. 

8. A method according to claim 7 wherein the measuring 
of the mean signal delays is accomplished in approximation 
by measuring phase delays of sinusoidally modulated sig- 
nals. 

9. A method according to claim 2 wherein said calculation 
includes the steps of 

i. determining matrix P given by 



P = 





2 


-Sn 


-Sn 




1 


2 


-s 2l 




-^23 


2 


2 




-s 32 


"^33 




2 


-S41 


-S42 


S 43 



the is the jth component of polarization launch S if 
ii. determining the inverse, P -1 , of matrix P, and 
hi. determining the four components of the vector 

t-(t 0 ,Q 1 ,Q 2 ,Q 3 ) where Q { is the ith component of the 

PMD vector Q from t^P -1 ^, where x g is a four- 
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component vector formed from the four measured of 
the mean signal delays. 

10. Amethod according to claim 9 wherein the measuring 
of the mean signal delays is accomplished in approximation 

5 by measuring phase delays of sinusoidally modulated sig- 
nals. 

11. Amethod according to claim 2 wherein the measuring 
of the mean signal delays is accomplished in approximation 

0 by measuring phase delays of sinusoidally modulated sig- 
nals. 

12. Amethod according to claim 1 wherein the measuring 
of the mean signal delays is accomplished in approximation 
by measuring phase delays of sinusoidally modulated sig- 

15 nals. 

13. A method according to claim 12 wherein the calcu- 
lation is based on the formula: 

2Q tanio^-To)— "qTtanCw^x/l) 

where co m is the angular modulation frequency, q is a unit 
vector pointing in the direction of the input polarization- 
mode-dispersion (PMD) vector, Ax is the magnitude of the 
25 PMD vector, x g is the mean signal delay as measured by the 
phase of a sinusoidally modulated signal, *c 0 is the intrinsic 
group delay of the device, and s is the input Stokes vector 
of a light signal transmitted through the device. 

14. A method according to claim 13 including 
30 determining, at said device output, the polarization state of 

each test light signal for determining the PMD of said device 
at said output thereof. 

15. A method according to claim 14 including using at 
35 least two of said test signals of known polarization state for 

determining the rotation matrix of said device. 

***** 
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